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1.1.  Introduction 

We  have  developed  a  model  for  the  generation  of  seismic  noise  at  the  sea  floor  which  will 
enable  us  to  extrapolate  measurements  made  at  the  sea  floor  with  ocean  bottom  seismographs  to 
depths  within  the  underlying  sediment  and  basement.  This  physical  model  is  consistent  with  sea 
floor  noise  measurements  we  have  made  at  Scripps  during  the  past  decade  and  largely  discounts 
the  hypothesis  that  noise  is  generated  at  the  sea  floor  by  current  motion  inducing  instrument 
’Vocking"  through  mechanisms  such  as  vortex  shedding.  We  have  never  been  able  to  demonstrate 
a  correlation  between  bottom  currents  and  seismic  noise  on  the  sea  floor  or  in  finely  controlled 
experiments  such  as  the  Lopez  Island  intercomparison  test.  However,  there  is  a  definite  correla¬ 
tion  between  sea  floor  noise  and  the  general  weather  in  the  vicinity  of  the  seismographs.  Details 
of  these  measurements  and  experiments  are  contained  in  previous  Annua)  and  Semi-Annual 
Technical  Reports. 

1.2.  Power  Spectra  of  Seafloor  Noise 

The  basic  assumptions  involved  in  the  model  are: 

(1)  Horizontal  stratification.  To  simplify  the  mathematics  involved  in  the  model  we  consider 
only  the  2-D  case,  where  field  variables  can  generally  depend  on  (r  ,z  ,t );  the  earth  structure 
depends  only  on  z  ,  however. 

(2)  Noise  sources  are  confined  to  the  water  column.  The  existing  data  are  consistent  with  the 
hypothesis  that  the  primary  sources  of  seismic  noise  on  the  sea  floor  are  hydrographic  distur¬ 
bances  in  the  water  column.  We  formalize  this  hypothesis  by  imposing  on  the  noise  model 
an  (unknown)  distribution  of  vertical  point  loads,  /  (r  ,t ),  at  the  sediment-water  interface 
(Figure  1).  Later,  in  the  fully  three-dimensional  calculations,  we  shall  let  /  depend  on  both 
r  and  6 ,  the  azimuth.  If  we  let  p(r,z,t)  be  the  incremental  pressure  field  in  the  water 
column  (z  <  0)  and  r(z,z,t )  be  the  stress  tensor,  then  the  (inhomogeneous)  boundary  con¬ 
ditions  to  be  applied  at  z  =  0  are  ru  (r  ,0,t )  =  /  (r,t)  =  -p(r,0,t)  and  r„  (r  ,0,1 )  =  0. 
These  boundary  conditions  have  the  desirable  effect  of  decoupling  the  problem  of  noise  in 
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the  ocean  from  noise  in  the  solid  medium  by  the  specification  of  /  (r  ,t ).  The  assumption 
that  noise  sources  are  confined  to  the  water  column  also  fixes  the  radiation  condition  for 
large  *  to  be  that  of  downgoing  waves  only.  In  our  model  this  condition  is  imposed  in  the 
half-space  underlying  the  layered  stack  representing  crustal  stratification  (Figure  1). 

(3)  The  noise  field  is  stationary  in  time  and  space  and  isotropic  in  horizontal  coordinates.  In 
the  2-D  problem  discussed  here  this  implies  that  its  autocorrelation  function  depends  only 
on  the  temporal  and  spatial  lag  variables: 

CH  (r.*)  =  <  /  ('.  ,t,  )f  ( |  r+r.  |  ,  t  +«.  )> 
oo  2r  oo 

=  f  dto  £r*drof  (ro>to).  /  (I  *+r.  I  »  ‘-I-*.)  (21) 

where  r  =  |  r|  ,  r0  =  |  r,  |  ,  r  r,  =  rrecos^„,  and  the  brackets  <  >  denote  an  ensemble 
average. 

These  assumptions  allow  a  precise  formulation  of  the  inverse  problem  for  C//  given  both 
seismic  and  acoustic  noise  observations  on  the  sea  floor;  they  also  allow  a  precise  formulation  of 
the  problem  of  predicting  the  noise  field  at  depth. 

Let  u(r  ,2  ,t )  be  the  vector-valued  displacement  field  at  some  position  z  >  0,  and  let 
G(|  x—r,  |,2,  t— 10)  be  the  displacement  response  at  (x  ,y  ,z  ,t)  to  a  vertical  point  load  at 
(*.  <Va  >0,*» ).  where  r2  =  z2+y2  and  r7  =  x7  +  y7.  Then  we  have 

U(r  ,e,t)  =  G(|  r-r0  |  ,  2 ,  t -t, )  *  /  (r  ,0  (2.2) 

where  the  asterisk  denotes  convolution  in  both  time  and  horizontal  space.  Because  the  source 
field  is  presumed  to  be  stationary  and  isotropic  and  G  depends  only  on  horizontal  and  temporal 
coordinate  differences,  the  tensor-valued  autorcorrelation  function  for  U  also  depends  only  on  the 
lag  variables: 

C„  (r  ,2  ,t )  =  <  U(rc  ,2  ,t0  )Ur(|  r-tr,  |  ,  2 ,  t  +t,  )>  (2.3) 

Since  /  ,  «  and  G  are  all  axisymmetric,  it  can  be  shown  (see  Appendix  Al.l)  that  equation  (2.3) 
is  equivalent  to: 
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c..(r.*,0  =  r(r,*,t)  *  C//  («•.*) 

where  F  is  the  tensor-valued  autocorrelation  function  for  G: 


oo  2  w  oo 

r(r  .*,*)  =  f  dl,  £r°  dr»G(r»>2>t)GrU  *0  +*1  ,  Z,  t„  +  t) 

If  the  Fourier-Besse)  transform  pair  A  (r  ,t )  and  h  ( k  ,w)  is  defined  by: 


00  00 

h(k,u)=  J  dt  t  f  rdrJ9  (hr  )h  (r  ,t ) 

-oo  *0 

00  00 

M'.O  -  -j t~f  d(Je~iU‘  fwfc/#(*r  )*(*,«) 

—  00  0 


(2.6b) 


then  the  transformation  of  C„  can  be  written  (see  Appendix  A1.3)  as: 


C..(t  ,2  ,W)  =  2*T(fc  l2  ,u>)C„  (k  ,u>)  (2.7) 

where  T(k  ,z  ,oi)  and  Ct  t  ( k  ,w)  denote  the  Fourier- Bessel  transforms  of  T(r  ,z  ,t )  and  C/ 1  (r  ,t ), 
respectively. 

A  special  case  of  C„  (r  ,u)  is  the  vertical  component  power  spectrum  P,  (cv),  which  is  the 
22 -component  of  C„  (r  ,ui)  at  r  =  0.  We  have: 

00 

P,(u>)  -  2z  £kdkr„(k,z,w)C/t  [k  ,u) 

00 

=  2x|w*|  G,(k,z,u)\7Cn(k,  u>),  (2.8) 

since  the  zz  component  of  T(Jt  ,z  ,u/)  is  G,*[k  ,w)G,  (k  ,w)  (see  Appendix  A1.3). 

The  spatial  variation  of  Ct /  is  poorly  constrained;  a  simply  parameterized  form  is 

C,,(r,U)  =  <?/,(«)«  (2  9) 

where  (to)  =  Cff  (o  ,to)  and  K*(w)  depend  on  frequency;  such  a  function  is  frequently  used  in 
studying  wave  propagation  through  inhomogeneous,  random  media  (Aki  and  Richards,  1980,  p. 

741).  The  Bessel  transform  of  e-**'*  is  — — j  e-**/4**,  so  that  for  such  a  C/t  ,  the  vertical 

2  K 


power  spectrum  is 


P,  («)  -  2*  C,a,  (to)  —1—  fkdk  |  G,  (k  ,z  ,to)  |  1  e  ->*/**'M  (2.10) 

2a  (u/J  *5 


G,(k,z  ,  u>)  can  be  computed  for  a  given  crustal  model  using  the  method  outlined  in  the  next  sec¬ 
tion  and  P,  (u>)  can  be  measured  at  the  sea  floor,  z  =  0;  Cj°j  (u>)  can  then  be  estimated  if  reason¬ 
able  values  of  K  (u>),  the  correlation  wavenumber,  arc  used.  Given  this,  the  noise  spectrum  at 
depth  z  >  0  can  be  predicted. 


1.3.  Green's  Functions 

The  viscoelastic  Green’s  functions  necessary  for  the  correlation  calculations  are  computed 
for  an  earth  model  consisting  of  a  stack  of  homogeneous  layers  including  a  fluid  ocean  as  the 
uppermost  layer.  The  problem  for  a  point  force  in  a  flat  earth  suggests  the  use  of  cylindrical  coor¬ 
dinates.  In  such  a  system,  the  displacements  can  be  expanded  in  terms  of  vector  cylindrical  har¬ 
monics: 


00  00 
V(r,8,z,t)  =  ~  f 

+  Ur(k,z,U)Sr+  Ut(k,z, u,)Tkm] 

where  r  ,  8  and  z  are  cylindrical  coordinates  and  R*",  S*m  and  T”  are  defined: 


Ykm  =  JM  (kr  )e  m  =  0,  ±  1,  ±  2, 
Rr(r  fi)  =  Y?e, 


T  r(r,6)  + 


1 

,  _L 

an"  . 

k 

TT  *' 

kr 

9  8  e* 

JL 

ZY”  . 

1 

9  Yf  . 

kr 

9  e  e' 

'  k 

9r  «* 

(3.1) 


(3.2) 


where  e, ,  t,  and  tk  are  unit  vectors  in  the  z ,  r  and  8  directions  (Takeuchi  and  Saito,  1972). 
For  a  vertical  point  force  the  problem  is  axisymmetric  and  m  =  0.  The  vertical  displacement  at 
r,  8  and  z  in  equation  (3.1)  reduces  to: 


V,  (r  ,8,z  ,w)  =  J^kdkU,  (k  ,z  ,u>)Ja  (kr  ) 


(3.3) 


In  performing  the  numerical  computations,  we  find  it  convenient  to  change  the  variable  of 
integration  to  horizontal  slowness  or  ray  parameter  p  ;  k  =  u/p  : 


The  physical  mode)  used  in  our  computations  is  shown  in  Figure  1.  In  this  case  an  ocean 
overlies  a  sequence  of  viscoelastic  layers  which  are,  in  turn,  underlain  by  a  homogeneous  half 
space.  Boundary  conditions  for  the  solution  include  a  traction-free  boundary  at  the  free  surface, 
continuity  of  traction  and  displacement  between  solid  layers,  continuity  of  ra  and  U,  at  the 
fluid-solid  boundary  and  a  requirement  that  waves  in  the  half  space  satisfy  the  radiation  condi¬ 
tions.  These  conditions  simply  require  that  only  downward  propagating  waves  exist  in  the  half¬ 
space.  The  vertical  force  at  the  sea  floor  can  be  thought  of  as  an  inhomogeneous  boundary  condi¬ 
tion  producing  a  jump  in  r„  at  the  ocean  bottom.  In  the  solid  layers  the  user  is  required  to 
specify  the  compressional  and  shear  velocities,  the  density  and  compressional  and  shear  Q ’s.  In 
the  fluid,  of  course,  only  the  acoustic  velocity,  density  and  compressional  Q  can  be  specified. 

We  compute  the  vertical  response  of  this  stack,  U,(p,z, tv),  with  a  technique  which  com¬ 
bines  phase-related  reflection  and  transmission  coefficients  in  an  iterative  fashion  (Apsel,  1978). 
The  method  is  essentially  identical  to  the  work  published  by  Kennett  and  Kerry  (1979)  and  Ken- 
nett  (1980).  The  fundamental  advantage  of  this  approach  is  the  achievement  of  numerical  stabil¬ 
ity  at  high  frequencies  and  wavenumbers. 

We  have  chosen  a  suite  of  four  simple  models  to  represent  the  structure  in  the  northwest 
Pacific  near  the  proposed  MSS  sites.  The  physical  properties  for  the  site,  see  Table  I,  were  taken 
from  the  Initial  Reports  of  the  Deep  Sea  Drilling  Project  for  Leg  32,  holes  303  and  304  in  the 
northwest  Pacific  (Figure  2).  The  models  encompass  the  range  of  sediment  thickness  and  types 
noted  in  the  Naval  Ocean  Research  and  Development  Activity  environmental  report  for  the 
Marine  Seismic  System  (MSS).  The  models  consist  of  approximately  5.5  km  of  ocean  overlying 
diatom  radiolarian  ooze  with  very  low  rigidity.  This  layer,  250-400  meters  in  thickness,  overlies 
either  a  hundred-meter  thick,  rigid  chert  layer  or  oceanic,  basaltic  basement. 

|  V,  [p  ,*,w)  |  for  a  400-meter  sedimentary  layer  overlying  basement  is  illustrated  in  Fig¬ 
ure  3  with  frequency  increasing  from  the  background  into  the  foreground  and  horizontal  slowness 
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increasing  from  left  to  right  across  the  foreground.  Figure  4  is  a  blowup  of  the  bracketed  section 
of  Figure  3.  The  dominant  pulse  in  the  frequency-slowness  spectra  corresponds  to  the  fundamen¬ 
tal  Stoneley  or  interface  wave  and,  at  high  frequency,  has  a  slowness  slightly  greater  than 
1/0.2  -5s  km-1.  At  low  frequencies  the  Stoneley  wave  begins  to  sense  the  underlying  chert  and 
basement  and  the  slowness  of  the  fundamental  decreases  (phase  velocity  increases);  that  is,  the 
mode  is  normally  dispersed.  We  note,  in  addition,  that  at  this  "knee"  below  1  Hz  in  the  funda¬ 
mental  dispersion  curve,  the  amplitudes  of  the  Stoneley  wave  decrease  rapidly  as  the  wave  enters 
increasingly  rigid  media.  At  small  slownesses  or  high  phase  velocities  small  ridges  correspond  to 
the  propagation  of  P  and  S  body  waves  and  leaky  modes.  In  Figure  4,  a  small  ridge  at  a  slow¬ 
ness  of  0.25  s  km-1  corresponds  to  the  compressional  head  wave  propagating  along  the  basement 
interface.  In  the  low  frequency  detail  plot  in  Figure  4,  we  see  the  fundamental  Stoneley  wave 
approaches  the  shear  wave  slowness  in  the  rock  (1/2.31  =  0.433  s  km-1). 

We  have  found  the  vertical  wavefunctions,  U2  (p  ,z  ,w),  for  various  propagation  modes  to  be 
very  useful  in  studying  the  effects  of  burial  of  the  MSS  seismometer.  A  model  for  thick  sediments 
including  a  chert  layer  overlying  basement  is  shown  in  Figure  5.  The  symbols  indicate  the 
receiver  burial  depths  at  which  we  shall  calculate  amplitude  spectra.  A  vertical  wavefunction, 
that  is  |  U,  (p  ,z  ,w)  |  as  a  function  of  depth,  for  the  fundamental  Stoneley  mode  branch  is 
shown  in  Figure  6  for  a  horizontal  slowness  of  4.5  s  km-1  and  a  frequency  of  0.2  Hz.  The  max¬ 
imum  displacement  occurs  at  the  sea  floor  where  a  discontinuity  marks  the  source  depth.  The 
wavefunction  decays  exponentially  into  the  overlying  water  and  into  the  underlying  chert  and 
basement.  The  mode  of  propagation  within  the  water  column  is,  of  course,  purely  acoustical.  For 
this  particular  important  mode  of  propagation,  the  effect  of  burial  is  clearly  to  rapidly  reduce 
vertical  displacement  and,  by  implication,  seismic  noise.  Figure  7  is  a  vertical  wavefunction  for 
the  first  overtone  of  the  fundamental  Stoneley  wave  at  a  frequency  of  1  Hz  and  a  horizontal  slow¬ 
ness  of  4  s  km-1.  In  this  case  there  are  two  oscillations  within  the  sediment  and  exponential 
decay  into  the  overlying  chert  and  basement.  In  this  case,  the  amplitude  of  this  mode  is  actually 
larger  toward  the  bottom  of  the  sediments,  but,  because  of  the  position  of  the  source,  the  mode  is 


not  a  significant  contributor  to  the  propagation  of  sea  floor  seismic  noise.  This  is  easily  seen  in 
Figure  3  where  the  fundamental  mode  clearly  dominates  any  of  the  overtones. 

We  have  studied  the  effect  of  burial  depth  at  all  frequencies  for  all  propagation  modes  by 
performing  the  zero-order  Hankel  Transform  in  equation  3.  The  complete  frequency-slowness 
spectrum  is  transformed  to  yield  frequency-distance  spectra  such  as  the  one  shown  in  Figure  8  for 
a  range  of  10  km.  The  squares  correspond  to  a  receiver  at  the  sea  floor  or  a  depth  of  5.4  km,  the 
triangles  are  200  meters  into  the  sediment,  the  diamonds  are  50  meters  into  the  chert,  the  +’s  50 
meters  into  the  basalt  and  finally  the  X’s  are  200  meters  into  the  basalt.  These  depths  are  shown 
graphically  on  the  model  in  Figure  5.  For  this  particular  case,  the  broad  peak  associated  with  the 
fundamental  Stoneley  mode  propagation  disappears  as  the  depth  of  burial  increases  from  within 
the  sedimentary  column  to  within  the  crystalline  basement.  Generally,  the  vertical  amplitudes 
decrease  monotonically  with  burial  depth  except  for  frequencies  in  excess  of  1-2  Hz  where  the  pat¬ 
tern  is  more  variable.  The  most  dramatic  effect  of  burial  depth  is  to  remove  the  broad  Stoneley 
peak  when  the  receivers  enter  the  more  rigid  chert  layer.  However,  we  are  reluctant  to  advocate 
burying  the  MSS  instrument  in  chert  as  a  solution  to  the  noise  problem.  We  have  modeled  the 
chert  as  a  homogeneous,  fairly  rigid  sedimentary  layer.  In  fact,  the  chert  exists  as  flat-lying, 
platy  deposits  within  a  less  rigid  sedimentary  matrix.  That  is,  the  chert  layer  is  most  properly 
modeled  as  a  composite  medium  in  which  the  rigidity  of  the  composite  may  vary  rapidly  in  both 
horizontal  and  vertical  directions.  Furthermore,  chert  is  extremely  difficult  to  drill  and  may 
prevent  penetration  to  depths  necessary  to  provide  'Shielding"  from  the  large  amplitude,  overlying 
sedimentary  Stoneley  waves. 

We  have  conducted  an  exhaustive  search  of  the  characteristics  of  wave  propagation  in  the 
four  representative  models  in  Table  I  with  a  vertical  point  force  at  the  sea  floor  as  a  source.  Fig¬ 
ures  9-35  are  amplitude  spectra  for  various  combinations  of  models  at  different  source-receiver 
separations  and  receiver  depths.  The  figures  are  quite  consistent  although  many  of  the  major  con¬ 
clusions  can  be  illustrated  by  an  examination  of  a  few  of  the  figures.  Figure  9  shows  amplitude 
spectra  for  thick  and  thin  sediments  (400  and  250  m  respectively)  each  with  a  hundred-meter 
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thick  chert  layer  below  the  sediment  and  above  the  basement.  We  have  chosen  a  range  of  10  km 
in  this  case  and  the  receivers  are  at  the  sea  floor  and  buried  50  meters  into  the  chert.  We  sec  the 
effect  of  burial  noted  before  in  the  discussion  of  Figure  8.  That  is,  burial  of  the  receiver  generally 
removes  the  broad  peak  between  0.1  and  2  Hz.  In  this  figure  we  note  the  peak  occurs  at  lower 
frequencies  and  is  greater  in  breadth  for  thick  sediments  as  compared  to  the  thin  sediments.  This 
is  certainly  reasonable  since  the  fundamental  Stoneley  mode  vertical  wavefunction  will  begin  to 
sense  the  rigid  chert  and  basement  at  higher  frequencies  in  the  case  of  the  thin  sediments. 

Figure  20  is  computed  for  the  same  pair  of  models  except  the  deep  receiver,  in  this  case,  is 
200  meters  into  the  basalt.  Again  the  broad  Stoneley  wave  peak  disappears  with  burial  in  the 
rigid  media  below  the  sediments.  We  note,  interestingly,  there  is  little  difference  in  the  spectra 
for  the  two  receivers  buried  within  the  basement.  We  might  expect,  therefore,  that  the  thickness 
of  the  unconsolidated  sediment  or  ooze  should  play  no  appreciable  role  in  the  selection  of  a  dril¬ 
ling  site  if  we  are  considering  the  seismic  noise  problem  alone.  The  would  not  be  true  of  an  OBS 
experiment  in  which  the  thinner  sediments  might  be  preferred  in  order  to  reduce  the  breadth  of 
the  Stoneley  peak. 

Figure  28  shows  the  effect  of  source-receiver  separation  on  the  point  force  amplitude  spectra. 
The  models  are,  again,  identical  to  those  already  discussed  and  the  receivers  are  at  the  sea  floor. 
As  the  range  increases  in  this  case  from  10  to  1100  km  the  Stoneley  peak  disappears  and  the 
overall  amplitude  drops  by  more  than  2  decades.  We  have  used  a  fairly  low  Q  for  the  sediments 
which  inevitably  leads  to  the  loss  of  the  sedimentary  Stoneley  waves  as  propagation  distance 
increases.  Geometrical  spreading  alone  would  decrease  amplitudes,  but  not  distort  the  shape  so 
severely  as  we  see  here.  Unfortunately,  we  know  very  little  of  the  shear  Q  in  sea  floor  sediments 
and  the  predictions  of  these  models  must  be  studied  with  this  difficulty  in  mind.  Certainly,  low 
group  velocity  (0.6  km  s~')  sedimentary  Stoneley  waves  have  been  observed  over  distances  of 
1000  km  on  the  sea  floor  (Okal  and  Talandier,  1980).  Further  research  on  this  phenomenon  is 
absolutely  necessary  and  will  be  considered  in  detail  in  our  future  work. 


Finally,  in  Figure  13  we  sec  the  effect  of  the  chert  layer  in  this  problem.  The  sedimentary 
layer  is  250  m  in  thickness  and  the  receivers  are  at  the  sea  floor  and  50  m  below  the  unconsoli¬ 
dated  sediments.  Very  clearly,  the  presence  or  absence  of  chert  in  the  bottom  of  the  sedimentary 


section  makes  little  difference  in  the  propagation  problem. 

1.4.  CONCLUSIONS 

We  have  completed  a  thorough  examination  of  the  Green’s  functions  used  to  propagate 
noise  in  our  ocean  bottom  seismic  noise  model.  Almost  certainly,  the  fundamental  model  Stone- 
ley  wave  is  responsible  for  transport  of  the  bulk  of  seismic  noise  energy.  We  can  make  a  number 
of  recommendations  at  this  point  with  respect  to  the  siting  of  the  MSS  system. 

1)  The  thickness  of  the  sedimentary  layer  is  not  important  in  site  selection  with  the  under¬ 
standing  the  instrument  will  be  installed  within  the  oceanic  basement. 

2)  The  presence  or  absence  of  a  thin  chert  layer  does  not  play  a  large  role  in  the  noise  levels 
we  expect  to  see  at  a  buried  receiver.  However,  as  noted  in  the  text,  even  though  the  fun¬ 
damental  Stoneley  wave  is  severely  attenuated  in  the  chert,  we  do  not  recommend  siting  at 
a  location  over  smooth  bottom  where  chert  is  likely  to  be  present.  The  drilling  ship,  cer¬ 
tainly,  will  have  an  unpredictable  amount  of  difficulty  in  penetrating  the  chert  layer. 

3)  Burial  of  the  seismic  sensor  in  the  basement  will  likely  reduce  seismic  noise  as  observed  at 
the  sea  floor  in  the  microseism  band,  0.1-2  Hz,  by  as  much  as  3  decades  in  displacement  and 
6  decades  (60  db)  in  power. 

4)  Because  surface  weather  controls  the  forcing  of  microseismic  noise,  the  sensor  should  be 
sited  in  an  area  with  weather  as  mild  as  possible. 

We  must  also  mention  some  precautions: 

1)  A  lack  of  understanding  of  the  shear  Q  could  affect  these  results  somewhat.  However, 
because  the  role  of  the  sedimentary  Stoneley  waves  is  fairly  unimportant  if  the  sensor  is 
buried  in  the  basement,  we  don’t  expect  the  effect  to  be  large. 


Table  1.  Model  Parameters 
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Figure  1.  An  illustration  depicting  some  of  the  assumptions  we  propose  to  use  in  formulating  a 
hypothetical  noise  model.  Noise  sources  within  the  water  column  are  represented  by  the 
vertical  load  distribution  they  impose  on  the  sea  floor,  /  (z,<).  In  the  half-space,  only 
downgoing  waves  are  allowed.  The  problem  is,  given  observations  of  seismic  and  acoustical 
noise  on  the  sea  floor  recorded  by  an  OHS,  to  predict  the  noise  spectrum  at  a  proposed  MSS 
site  at  depth. 
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Figure  2.  Chart  of  northwest  Pacific  and  contiguous  masses  showing  proposed  Sites  1  and  2  for 
MSS  emplacement.  Sedimentary  thicknesses  arc  greater  at  Site  2  although  the  presence  of 
chert  is  not  as  ubiquitous  as  at  Site  I.  The  weather  is  generally  better  at  Site  2  although 
wave  heights,  an  important  consideration  for  seismic  noise,  are  generally  quite  large  year 
around  in  this  part  of  the  Pacific. 
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Figure  3.  Frequency  vs.  horizontal  slowness  for  vertical  displacement  at  the  sea  floor  due  to  a 
vertical  point  force  at  that  depth.  The  function  depicted  is,  in  fact,  the  modulus  of  the 
complex  vertical  displacement. 


TRANSFORMED  VERTICAL  DISPLACEMENT  DUE  TO  A  VERTICAL  POINT  FORCE 


SOURCE  AND  RECEIVER  AT  THE  SEDIMENT-WATER  INTERFACE 


Figure  4.  Blowup  of  low  frequency,  low  horizontal  slowness  for  the  modulus  of  vertical  displace¬ 
ment  at  the  sea  floor  due  to  a  vertical  point  force  at  that  depth.  The  area  of  this  plot 
corresponds  to  the  area  of  the  brackets  in  Figure  S. 
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Figure  5.  Compression&l  (dashed  line)  and  shear  velocity  (solid  line)  structure  for  a  typical  MSS 
site.  400  meters  of  ooze  with  very  low  rigidity  (shear  velocity  =  200  ms-1)  overlies 
100  meters  of  fairly  competent,  high-velocity  chert.  The  chert,  finally,  overlies  ocean  base¬ 
ment  rocks  with  a  compressional  velocity  of  4  km  s'1.  The  density  of  ooze  is  only  slightly 
greater  than  that  of  the  overlying  ocean.  The  arrows  and  symbols  denote  receiver  depths 
for  which  vertical  displacement  amplitude  spectra  were  computed  (sec  Figure  8)  for  a  source 
at  the  sea  floor. 
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Figure  6.  |  U,  ( p  ,w,z  )  |  plotted  versus  depth  for  a  frequency  of  0.2  Hz  and  a  horizontal  slowness 
of  4.5  s  km'1.  This  point  in  Figure  3  corresponds  to  the  fundamental  Stonelcy  mode.  The 
discontinuity  at  a  depth  of  5.4  km  is  due  to  the  presence  of  the  source  at  the  sea  floor. 
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Figure  8.  Spectra  at  a  source-receiver  separation  of  10  km  for  all  modes  of  propagation.  The 
receiver  depths  are  keyed  to  the  arrows  and  symbols  in  Figure  5.  The  two  deepest  receivers 
are  within  the  oceanic  basement  rocks.  The  broad  peaks  for  the  sea  floor  and  mid¬ 
sedimentary  receivers  arc  caused  almost  solely  by  the  propagation  of  low  phase  and  group 
velocity  Stoneley  waves.  We  note  the  dramatic  decrease  in  amplitude  as  the  receiver  moves 
below  the  sedimentary  ooze  layer  into  the  chert  and  underlying  basement.  The  spectra  for 
the  shallow  receivers  is  obviously  very  band-limited  and,  although  not  shown  here,  the  peak 
frequency  increases  as  the  sedimentary  thickness  decreases.  The  spectra  were  obtained  by 
computing  a  zero-order  liankcl  transform  of  several  depth  points  on  vertical  wavcfunctions 
such  as  that  in  Figure  C. 
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Figure  9.  Amplitude  spectra  at  10  km  range  for  two  crustal  models  consisting  of  a  100-meter 
thick  chert  layer  over  a  basement  and  below  either  a  thick  (400  m)  or  a  thin  (2.10  in)  layer 
of  sediments;  overlying  the  sediments  is  5.4  km  of  fluid.  The  source  is  a  vertical  point  force 
at  the  sea  floor  and  receivers  are  at  the  sea  floor  and  50  meters  into  the  chert. 
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Hgure  10.  Amplitude  spectra  at.  1100  km  range  for  two  crustal  models  consisting  of  a  100-m 
thick  chert  layer  over  a  basement  and  below  either  a  thick  (400  m)  or  a  thin  (250  in)  layer 
of  sediments;  overlying  the  sediments  is  5.4  km  of  fluid.  The  source  is  a  vertical  point  force 
at  the  sea  floor  and  receivers  are  at  the  sea  floor  and  50  meters  into  the  chert. 
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Figure  13.  Amplitude  spectra  at  10  km  range  for  two  crustal  models  consisting  of  a  thin  (250  m) 
sediment  layer  overlain  by  5.4  km  of  fluid  and  above  either  a  basement  or  100  meters  of 
chert  and  a  basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and  receivers  arc 
at  the  sea  floor  and  50  meters  below  the  sediments. 
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Figure  14.  Amplitude  spectra  at  1100  km  range  for  two  crustal  models  consisting  of  a  thin 
(250  m)  sediment  layer  overlain  by  5.4  km  of  fluid  and  above  cither  a  basement  or 
100  meters  of  chert  and  a  basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and 
receivers  are  at  the  sea  floor  and  50  meters  below  the  sediments. 
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Figure  15.  Amplitude  spectm  at  10  km  range  for  two  crustal  models  consisting  of  a  thick  (400  m) 
sediment  layer  overlain  by  5.4  km  of  fluid  and  above  cither  a  basement  or  100  meters  of 
chert  and  a  basement.  1  he  source  is  a  vertical  point  force  at  the  sea  floor  and  receivers  are 
at  the  sea  floor  and  50  meters  below  the  sediments. 
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ure  16.  Amplitude  spectra  at  1100  km  range  for  two  crustal  models  consisting  of  a  thick 
(400  m)  sediment  layer  overlain  by  5.4  km  of  fluid  and  above  either  a  basement  or 
100  meters  of  chert  and  a  basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and 
receivers  arc  at  the  sea  floor  and  50  meters  below  the  sediments. 
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Figure  18.  Amplitude  spectra  at  ranges  10  and  1100  km  for  two  crustal  models  consisting  of  a 
thick  (400  m)  layer  of  sediments  overlain  by  5.4  km  of  fluid  and  above  either  a  basement  or 
100  meters  of  chert  and  a  basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and 
the  receivers  are  50  meters  below  the  sediments. 
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Hguro  19.  Amplitude  spectra  at  ranges  10  and  1100  km  Tor  two  crustal  models  consisting  of 
either  a  thin  (250  m)  or  thick  (400  rn)  sediment  layer  above  a  basement  and  below  5.4  km 
of  fluid.  1  he  source  is  a  vertical  point  force  at  the  sea  floor  and  the  receivers  are  50  meters 
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Figure  23.  Amplitude  spectra  at  1100  km  range  for  two  crustal  models  consisting  of  eilher  a  thick 
(400  m)  or  thin  (250  m)  sediment  layer  above  a  more  rigid  basement  and  below  5.4  km  of 
fluid.  The  source  is  a  vertical  point  force  at  the  sea  floor  and  receivers  are  at  the  sea  floor 
and  200  meters  into  the  basement. 
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Figure  24.  Amplitude  spectra  at  10  km  range  for  two  crustal  models  consisting  of  a  thin  (250  ni) 
sediment  layer  overlain  t>y  5.4  km  of  fluid  and  aln>ve  either  a  basement  or  100  meters  of 
chert  and  a  basement  The  source  is  a  vertical  point  force  at  the  sea  floor  and  receivers  are 
at  the  sea  floor  and  300  meters  below  the  sediments,  in  the  basement. 
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Figure  25.  Amplitude  spectra  at  1100  km  range  for  two  crustal  models  consisting  of  a  thin 
(250  m)  sediment  layer  overlain  by  5.4  km  of  fluid  and  above  either  a  basement  or 
100  meters  of  chert  and  a  basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and 
receivers  arc  at  the  sea  floor  and  300  meters  below  the  sediments,  in  the  basement. 
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Figure  2G.  Amplitude  spectra  at  10  km  range  for  two  crustal  models  consisting  of  a  thick  (400  m) 
sediment  layer  overlain  by  5.4  km  of  fluid  and  above  either  a  basement  or  100  meters  of 
chert  and  a  basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and  receivers  arc 
at  the  sea  floor  and  300  meters  below  the  sediments,  in  the  basement. 
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Figure  27.  Amplitude  spectra  at  1100  km  range  for  two  crustal  models  consisting  of  a  thick 
(400  m)  sediment  layer  overlain  by  5.4  km  of  fluid  and  above  cither  a  basement  or  100 
meters  of  chert  and  a  basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and 
receivers  are  at  the  sea  floor  and  300  meters  below  the  sediments,  in  the  basement. 
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f  igure  30.  Amplitude  spectra  at  ranges  10  and  1100  km  for  two  crustal  models  consisting  of 
either  a  thin  (250  m)  or  thick  (400  rn)  sediment  layer  overlain  by  5.4  km  of  fluid  and  a 
basement.  The  source  is  a  vertical  point  force  and  the  receivers  are  at  the  sea  floor. 
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Figure  31.  Amplitude  spectra  at  ranges  10  and  1100  km  for  two  crustsl  models  consisting  of 
either  a  thin  (250  m)  or  thick  (400  in)  sediment  layer  overlain  by  5.4  km  of  fluid  and  a 
basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and  the  receivers  are  200 
meters  into  the  basement. 
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Figure  33.  Amplitude  spectra  at  ranges  10  and  1100  km  for  two  crustal  models  consisting  of  a 
thin  (250  m)  layer  of  sediments  overlain  by  5.4  km  of  fluid  and  above  cither  a  basement  or 
100  meters  of  chert  and  a  basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and 
the  receivers  are  300  meters  below  the  sediments,  in  the  basement. 
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Figure  34.  Amplitude  spectra  at  ranges  10  and  1100  km  for  two  crustal  models  consisting  of  a 
thick  (400  m)  layer  of  sediments  overlain  by  5.4  km  of  fluid  and  above  either  a  basement  or 
100  meters  of  chert  and  a  basement.  The  source,  a  vertical  point  force,  and  the  receivers 
arc  at  the  sea  floor. 
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Figure  35.  Amplitude  spectra  at  ranges  10  and  1100  km  for  two  crustal  models  consisting  of  a 
thick  (400  mj  layer  of  sediments  overlain  by  5.4  km  of  fluid  and  above  either  a  basement  or 
100  meters  of  chert  and  a  basement.  The  source  is  a  vertical  point  force  at  the  sea  floor  and 
the  receivers  are  300  meters  below  the  sediments,  in  the  basement. 
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CHAPTER  II 


ANALYSIS  OF  AMBIENT  SEISMIC  NOISE  RECORDED  BY 
DOWNHOLE  AND  OCEAN  BOTTOM  SEISMOMETERS  ON 
DEEP  SEA  DRILLING  PROJECT  LEG  <8B 
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RICHARD  G.  ADAIR 
JOHN  A.  ORCUTT 
THOMAS  H.  JORDAN 


II. 1.  Introduction 


The  measurement  of  ambient  seismic  noise  near  the  sea  floor  is  requisite  to  formulating 
optimal  strategies  for  the  deployment  of  marine  seismometer  systems.  There  are,  however,  no 
published  studies  of  noise  within  the  oceanic  crust,  and  few  of  noise  at  the  sea  floor,  especially  of 
reliable,  absolute  noise  levels.  Most  available  measurements  were  made  during  the  VELA  Uni¬ 
form  Project  sponsored  by  the  Advanced  Research  Projects  Agency  in  the  1960’s.  It  was  hoped 
that  noise  levels  at  the  sea  floor  would  be  comparable  to  those  on  land  (Prentiss  and  Ewing, 
1963),  but  in  fact  they  were  found  to  be  significantly  higher  (Bradner  and  Dodds,  1964;  Schneider 
and  Backus,  1964;  Schneider  et  al- ,  1964;  Latham  and  Sutton,  1966;  Latham  and  Nowroozi,  1968). 
These  studies  suggest  that  the  noise  propagates  as  an  evanescent  wave  (Stoneley  wave)  trapped 
near  the  ocean  bottom-sea  floor  interface  (Bradner  et  al.,  1965;  Latham  and  Sutton,  1966;  Latham 
and  Nowroozi,  1968),  and  thus  reduced  noise  levels  are  expected  within  the  oceanic  crust.  The 
technology  and  expertise  needed  to  implant  instruments  in  the  sea  floor  has  only  recently 
developed,  primarily  through  the  Deep  Sea  Drilling  Project  (DSDP).  Borehole  seismometers  have 
been  used  in  oblique  seismic  experiments  on  several  DSDP  legs  (Stephen,  1980;  Stephen  et  al., 
1981;  Stephen,  1982),  but  the  data  are  not  suited  for  the  quantitative  analysis  of  seismic  noise.  A 
downhole  seismometer  system  sponsored  by  the  Defense  Advanced  Research  Projects  Agency,  the 
Marine  Seismic  System  (MSS),  was  deployed  in  Hole  395A  (Figure  1)  during  Leg  78B  with  the 
intent  of  recording  high-quality  noise  and  seismic  data.  This  paper  presents  the  first  quantitative 
characterization  of  ambient  seismic  noise  within  the  oceanic  crust.  In  addition,  a  direct  com¬ 
parison  of  noise  characteristics  at  and  below  the  sea  floor  is  made  using  data  collected  on  a  proxi¬ 
mate  ocean  bottom  seismograph  (OBS)  (see  Figure  1). 

II. 2.  Instrumentation  and  Operations 

The  borehole  instrumentation  of  the  MSS  consisted  of  state-of-health  sensors  and  two 
vertical-component,  short-period  seismometers  (Teledyne  Geotech  model  S-700'),  one  configured 
approximately  50  cm  below  the  other.  Each  seismometer’s  output  was  split  into  three  channels 
and  the  gains  of  the  resultant  six  channels  were  staggered  with  overlap  to  yield  a  dynamic  range 


.v 


■c. 


n..fcW*«.W 


jLXm 


of  144  dB  (Figure  2).  The  data  streams  were  filtered  to  prevent  aliasing,  digitized  at  75  samples 
per  second  using  a  10  bit  digital  word,  and  transmitted  via  cable  to  a  shipboard  recorder.  The 
cable  was  also  used  to  recover  the  borehole  instrumentation  package.  The  computed  displace¬ 
ment  response  for  the  combined  seismometer  and  filter  stages  is  shown  in  Figure  3.  The  response 
is  peaked  near  11  Hz,  and  rolls  off  at  12  dB/octave  towards  lower  frequencies  and  at  15 
dB/octave  to  the  Nyquist  frequency. 

The  borehole  package  was  lowered  from  Glomar  Challenger  to  a  depth  of  609  meters  in  Hole 
395A  (517  meters  within  the  basement),  where  cave-ins  blocked  access  to  the  bottom  42  meters. 
The  package  was  in  place  for  29.5  hours  but  data  were  recorded  only  during  the  final  26  hours, 
which  included  a  refraction  experiment  (see  Figure  4).  Although  the  seismometers  were  not 
clamped  in  the  hole,  the  adjacent  unfractured  and  massive  wall  rock  doubtlessly  promoted  proper 
coupling.  (The  wall  rock  character  was  inferred  from  caliper  and  acoustic  reflection  logs  con¬ 
ducted  earlier  during  Leg  78B.)  Ship  motions  conveyed  to  the  package  were  quelled  with  0.7  km 
of  slack  cable  payed  out  on  the  sea  floor.  Noise  samples  were  recorded  between  refraction  shots 
and  during  breaks  in  the  refraction  experiment,  both  intentional  and  unintentional  (when  explo¬ 
sives  failed  to  detonate).  The  times  of  relevant  operations  and  events  are  charted  in  Figure  4. 

Four  triaxial  ocean  bottom  seismographs  (OBS’s)  were  dropped  from  USNS  Lynch  by  the 
group  from  the  University  of  Texas  (UT)  Institute  of  Geophysics.  Each  OBS  was  programmed  to 
record  40  seconds  of  digital  data  with  a  sampling  interval  of  8.352  milliseconds  during  refraction 
shots  and  the  intentional  shooting  hiatus.  The  OBS  seismometers  were  10  Hz  Mark  L6  geo¬ 
phones.  1  he  OBS  displacement  response  (Figure  3)  increases  with  frequency  at  18  dB/octave  up 
to  10  Hz,  and  6  dB/octave  thereafter.  A  variety  of  malfunctions  prevented  data  recovery  except 
from  the  vertical  component  of  an  OBS  located  approximately  1-2  km  west  of  the  borehole  (see 
figure  1  for  launch  point).  A  more  precise  location  based  on  water  wave  arrivals  is  not  possible 
because  recordings  of  explosion  signals  but  not  of  noise  were  severely  clipped  by  frame  resonances 
above  10  Hz.  Eleven  34-second  samples  of  noise  were  obtained  from  the  OBS  data. 
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11.3.  Data  Quality 

Intrinsic  earth  noise  was  distinguished  from  instrument  and  ship  noise  on  the  basis  of  tem¬ 
poral  and  spectral  characteristics.  Spectral  characteristics  were  determined  from  power  densities 
and  coherences  computed  in  a  manner  described  by  Welch  (1967).  Time  series  were  broken  into 
segments  which  were  ‘de-meaned’,  tapered  with  a  normalized  Hanning  window,  and  Fourier 
transformed.  Auto-  (power)  and  cross-spectra  were  averaged  (stacked)  to  reduce  spectral  estimate 
variance.  Stable  estimates  were  generally  achieved  with  stacks  of  eight  segments,  although  as 
many  as  twenty-four  segments  were  used  when  available. 

Of  the  six  MSS  data  channels  (see  Figure  2),  only  three  were  usable  for  noise  analysis  (Ll, 
U2,  L2);  the  highest-gain  channel  (Ul)  was  always  clipped,  and  the  two  of  lowest  gain  (U3,  L3) 
were  dominated  by  least-count  noise.  An  example  of  MSS  data  from  Channel  Ll  is  shown  in  Fig¬ 
ure  5.  Typically,  the  data  were  dominated  by  monotonous,  1-4  second  energy  while  equally  per¬ 
sistent  high  frequency  energy  (>  10  Hz)  was  of  much  lower  amplitude.  Features  in  Figure  5 

which  are  attributed  to  instrument  noise  or  not  considered  to  represent  ground  motion  are  the 
dropped  bit  near  11.5  seconds;  the  frequent,  random,  unidirectional  (sensor  mass  up)  transients  of 
variable  magnitude  (e.g.  near  10.5  seconds);  and  the  truncated,  low-frequency  peak  centered  at  6 
seconds  (a  rare  event).  Single,  dropped  bits,  clearly  due  to  instrument  errors,  were  replaced  with 
linearly-interpolated  values.  High-frequency  amplitudes  and  occurrence  rates  of  the  transients 
generally  decreased  with  time,  while  the  low-frequency  amplitudes  increased.  Power  densities  of 
MSS  data  and  instrument  noise  are  shown  in  Figure  6.  The  data  were  eight  13.65  second  samples 
from  Channel  Ll  not  corrected  for  instrument  response.  Electronics  noise  levels,  available  only 
for  frequencies  greater  than  0.3  Hz,  were  estimated  from  RMS  noise  power  values  of  one-half 
octave  bandwidth  provided  by  Teledyne  Gcotech.  The  least-count  noise  level  is  that  expected 
from  uniformly-distributed  errors  due  to  the  representation  of  a  continuous  variable  by  a  finite, 
digital  word  (e.g.,  Otnes  and  Enochson,  1972).  Neither  of  these  sources  of  instrument  noise  dom¬ 
inated  above  0.3  Hz.  Least-count  noise  became  increasingly  important  at  frequencies  where  the 
instrument  response  rolled  off  and  resulted  in  diminished  coherence  of  energy  between  the  two 
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MSS  seismometers.  Power  densities  at  frequencies  below  O.lti  1 1  z  are  inferred  to  be  dominated  by 
instrument  noise  on  the  basis  of  low  coherence  (Figure  8) 

The  transients  seen  in  Figure  5  are  assumed  to  be  due  either  to  movement  of  the  unclamped 
borehole  package  or  to  pressure  impulses  caused  by  the  escape  of  the  casing  gas  (helium).  The 
transients  strongly  resemble  the  MSS  response  to  a  step  in  displacement  (Figure  7)  and  are  taken 
as  such  for  the  purpose  of  determining  power  densities.  The  transients  were  modeled  as  a  com¬ 
pound  Poisson  process  (Appendix  II),  and  their  predicted  power  density  P(j  )  at  frequency  /  is 

P  (}  )  ~  2A<  o2>  |  tf(/)|  2 

where  A  is  the  mean  rate  of  transient  occurrence,  <  a2>  is  the  mean  squared  amplitude,  and 
«{!)  is  the  Fourier  transform  of  a  transient  of  unit  amplitude.  A  and  <  a2>  were  estimated  for 
the  data  whose  power  densities  are  shown  in  Figure  6  and  the  resulting  P(f)  is  plotted  in  the  same 
figure  as  the  curve  labeled  ‘Transient  Noise’.  The  remarkably  close  agreement  between  the 
observed  and  predicted  noise  levels  implies  that  energy  beyond  the  break  in  power  density  slope 
(near  2.2  Hz)  is  concentrated  almost  entirely  in  the  transients. 

System  noise  levels  are  not  known  for  the  OBS,  nor  can  coherences  be  estimated  since  only 
the  vertical  component  of  one  instrument  was  functional.  Frequencies  at  which  instrument  noise 
dominated  were  qualitatively  determined  from  uncorrected  power  densities  (Figure  9).  The  OBS 
frame  resonance  is  dominant  between  10  and  25  Hz,  and  the  essentially  flat  portions  of  the  spec¬ 
trum  beyond  4  Hz  and  below  0.18  Hz  imply  instrumental  sources  of  white  noise.  In  addition,  esti¬ 
mate  stability  could  not  be  achieved  below  0.18  Hz  with  the  available  time  series  lengths. 

Ship-generated  acoustic  signals  undoubtedly  contributed  to  local  microscismic  noise. 
Sources  aboard  Challenger  include  machinery  and  the  dynamic  positioning  system.  The  dynamic 
positioning  system  was  silenced  only  once  during  the  MSS  deployment  (0101-01 20Z,  29  March), 
but  as  this  was  prior  to  the  start  of  data  recording,  a  quantitative  assessment  of  its  effects  is  not 
possible.  Unpublished  ocean  bottom  noise  spectra  estimated  from  data  collected  on  Scripps  Insti¬ 
tution  of  Oceanography  OBS’s  suggest  that  noise  caused  by  nearby  ships  is  significant  only  above 


-  16  - 


approximately  5  Hz  and  is  typically  manifested  as  sharp  peaks  at  a  few  frequencies  and  several  of 
their  higher  harmonics.  The  high  noise  levels  above  10  Hz  on  the  UT  OBS  may  be,  at  least  in 
part,  attributed  to  ship  noise  amplified  by  frame  resonances. 

Another  possible  noise  source  of  is  the  flow  of  water  past  the  instruments  (Duennebier  et  al., 
1981).  This  noise  is  thought  to  be  generated  when  current  speeds  are  in  excess  of  10  cm/scc,  but 
there  is  a  complex  dependence  on  body  geometry.  In  the  simple  case  of  a  cylinder,  the  frequency 
of  the  noise  is  inversely  proportional  to  the  diameter;  the  whip  antenna  on  the  OBS  would  have  a 
flow-noise  frequency  of  about  6  Hz  for  a  current  speed  of  10  cm/sec.  Near-bottom  current  meas¬ 
urements  (5  meters  above  the  sea  floor)  were  made  approximately  4  km  north  of  the  UT  OBS 
during  the  period  28  March-1  April  (Boyd,  this  volume).  The  average  current  speed  was  4.4 
cm/sec  and  the  maximum  was  8.9  cm/sec,  so  it  is  likely  that  current  noise  is  not  significant  for 
the  OBS  noise  observations.  In  the  case  of  the  MSS,  temperature  observations  indicate  that 
water  flowed  down  the  hole  at  the  rate  of  25  m/hr  and  had  disappeared  into  the  basement  within 
500  meters  sub-bottom,  100  meters  above  the  borehole  instrument  package.  Also,  there  was  no 
evidence  for  hydrothermal  circulation  anywhere  in  the  hole.  Hence,  water  flow  is  not  thought  to 
be  a  large  contributor  to  MSS  noise  measurements.  Thus,  we  feel  the  MSS  data  reliably  represent 
intrinsic  earth  noise  in  the  band  0.16-2.2  Hz,  while  the  OBS  data  probably  do  so  between  0.18 
and  4  Hz. 

II.4.  Discussion 

While  stability  of  noise  estimates  computed  from  44  minutes  of  MSS  data  implies  quasi- 
stationarity  on  the  time  scale  of  one  hour,  recorded  amplitudes  increased  enough  over  one  day  to 
cause  clipping  on  the  channel  of  highest  usable  gain  (Ll).  On  the  time  scale  of  ten  hours,  power 
densities,  computed  from  Channel  U2  data,  were  stable  in  shape  but  increased  by  up  to  5  dB 
(Figure  10).  Displacement  spectra  were  peaked  between  0.2  and  0.S  Hz  and  decreased  with  fre¬ 
quency  at  80  dB/dccade.  The  earliest  of  the  spectra  in  Figure  10  had  a  peak  value  of  4  x  10® 
nmJ/Hz  at  0.21  Hz  and  fell  from  1  x  10®  nmJ/Hz  at  0.33  Hz  to  1  nmJ/Hz  at  1.9  Hz.  A  more 
detailed  study  of  the  temporal  noise  variation  was  conducted  using  integrated,  uncorrcctcd  power 
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densities  from  Channel  U2  (which  was  not  clipped  by  noise)  as  a  measure  of  level  (Figure  lib). 
The  power  densities  were  computed  from  13.65  seconds  of  data  sampled  at  approximately  15 
minute  intervals;  these  realizations  of  power  were  smoothed  with  a  three-point  running  mean.  A 
trend  of  increasing  noise  levels  is  clearly  discernible  in  Figure  lib.  This  variability  of  ocean  bot¬ 
tom  noise  may  be  caused  by  changes  in  local  sea  state  and/or  weather  (Latham  and  Sutton,  1966; 
Latham  and  Nowroozi,  1968).  Using  bi-hourly  observations  recorded  in  Challenger’s  deck  log, 
only  swell  height  (Figure  11a)  was  strongly  correlated  with  noise  level.  The  (nonlinear)  interaction 
between  opposing  swell  produces  significant  pressure  fluctuations  on  the  sea  floor  at  half  the  swell 
period  (Longuet-Higgins,  1950;  Kadota  and  Labianca,  1981),  and,  indeed,  the  reported  swell 
period  of  10  -  11  seconds  is  approximately  twice  that  of  the  spectral  peak  between  3  and  5  seconds 
(Figure  10).  However,  a  more  definite  conclusion  is  not  justified  by  the  limited  data  set  use  here. 
The  period  of  the  noise  peak  is  within  the  range  of  the  so-called  "microseism  peak"  of  noise 
observed  at  the  ocean  bottom  (Latham  and  Sutton,  1966;  Latham  and  Nowroozi,  1968)  and  on 
continents  (Melton,  1976),  where,  in  the  latter  instance,  there  is  ample  evidence  ascribing  the 
microseism  peak  to  the  nonlinear  swell  interaction  (Haubrich  et  al.,  1963;  Haubrich  and  McCamy, 
1969). 

The  detailed,  temporal  behavior  of  noise  at  the  sea  floor  during  the  MSS  deployment  cannot 
be  reliably  determined  because  of  the  paucity  of  appropriate  OBS  data.  However,  simultaneous 
noise  observations  were  made  by  the  two  seismic  systems  on  eleven  occasions,  nine  of  which  may 
be  used  for  comparison  (Figure  4).  Six  of  these  were  recorded  over  a  75-minute  period  scheduled 
specifically  for  noise  study  and  afford  stable,  high-resolution  spectral  estimates  (Figure  12).  In  the 
band  of  overlap  (0.18-2  Hz),  noise  levels  were  higher  at  the  sea  floor  than  in  the  borehole  by  an 
amount  which  increased  with  frequency  from  10  dB  at  0.2  Hz  to  28  dB  at  2.2  Hz.  This  is  con¬ 
sistent  with  the  notion  that,  at  these  frequencies,  microseismic  noise  propagates  primarily  as 
fundamental-mode  Stoneley  waves  (Latham  and  Sutton,  1966;  Latham  and  Nowroozi,  1968; 
Brandner  and  Dobbs,  1964;  Bradncr  et  al.,  1965;  Tuthill  et  al.,  1981).  A  large  acoustic-shear 
impedance  mismatch  at  a  fluid-solid  interface  results  in  the  severe  exponential  decay  of  this  wave 
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away  from  the  interface,  essentially  trapping  the  seismic  energy  (Latham  and  Sutton,  1966; 
Tuthil!  et  al.,  1981).  Site  395  is  a  ponded  depression  in  the  midst  of  rugged  bathymetric  relief 
(Figure  1),  which  may  affect  seismic  amplitudes  in  complex  ways,  but  if  the  sediments  are  the 
dominating  influence  at  this  site,  Stoneley-wave  amplitudes  at  frequencies  above  that  of  the 
microseismic  peak  should  be  significantly  lower  below  the  base  of  the  sediments  than  at  the  sea 
floor. 

The  remaining  five  noise  windows  were  widely  spaced  in  time  and  were  intended  to  be  for 
explosion  signals,  but  the  shots  failed  to  detonate;  only  three  may  be  used  for  comparison  because 
of  excessive  instrument  noise  on  the  MSS.  Sea  floor  noise  levels  were  10-30  dB  higher  than  in  the 
borehole  for  these  samples,  as  well.  (The  power  densities  of  Figure  13  were  of  similar  resolution 
and  variance:  the  MSS  data  were  eight  6.83-second  samples  and  the  UT  OBS  data  were  eight 
4.28-second  samples.) 

Representative  extrema  of  microseismic  noise  levels  observed  during  the  ten  years  of  the 
Scripps  Institution  of  Oceanography’s  OBS  program  are  compared  in  Figure  14  with  the  lowest, 
reliable  levels  observed  during  Leg  78B.  The  SIO  OBS’s  employed  1-Hz  digitally  recorded 
seismometers  (Prothero,  1974;  Moore  et  al.,  1981).  Power  densities  were  computed  from  the 
vertical-component  output  which  resided  in  semiconductor  memory  prior  to  being  buffered  to  tape 
and  thus  were  not  contaminated  by  the  mechanical  tape-recorder  noise.  Data  for  the  ‘noisiest  site’ 
were  fifty-one,  8-second  samples  collected  southeast  of  Hawaii,  and  those  for  the  ‘quietest  site’ 
were  thirty-four,  4-second  samples  collected  on  the  East  Pacifc  Rise  near  20"  N.  The  extrema  are 
separated  by  20-30  dB.  The  UT  OBS  values  are  similar  to  the  SIO  OBS  maxima  while  the  MSS 
values  are  intermediate  to  low. 

Continental  noise  observations  are  compared  with  the  MSS  spectrum  in  Figure  15.  The 
Queen  Creek  data  were  recorded  in  a  mine  in  Queen  Creek,  Arizona  (Melton,  1976),  and  the  Laji- 
tas  data  were  recorded  at  a  site  near  Lajitas,  Texas  (Herrin,  1982).  The  Brune-Oliver  curves 
represent  average  and  high  levels  for  continental  stations,  and  were  taken  from  Melton  (1976), 
who  reanalyzed  Brune  and  Oliver’s  (1959)  summary  of  published  noise  observations.  MSS  levels 
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are  35-40  dB  greater  than  the  exceptionally  low  Queen  Creek  and  Lajitas  values,  and  fall  between 
the  Brune-Oliver  curves. 

II. 5.  Conclusions 

The  noise  levels  observed  by  the  MSS  in  Hole  595A,  approximately  500  m  below  the 
sediment-water  interface,  were  10-30  dB  lower  than  those  recorded  by  the  University  of  Texas 
OBS  over  the  frequency  range  0. 2-2.0  Hz  (Figures  12  and  13).  This  observation  is  consistent  with 
the  hypothesis  that  seismic  noise  at  these  frequencies  propagates  primarily  as  low-order  Stoneley 
modes  trapped  at  the  sea  floor  interface  (Bradner  et  al.,  1965;  Latham  and  Sutton,  1966;  Latham 
and  Nowroozi,  1968;  Tuthill  et  al.,  1981).  The  correlation  of  the  noise  levels  with  swell  heights 
(Figure  11)  suggests  that  the  generation  of  sea-floor  noise  at  this  mid-Atlantic  site  is  dominated 
by  local  nonlinear  ocean  wave  interactions.  To  within  the  accuracy  of  the  deck-log  observations, 
the  swell  period  was  approximately  twice  that  of  the  microseismic  noise  peak,  consistent  with  the 
Longuet-Higgins  mechanism  and  previous  ocean  bottom  measurements  (Latham  and  Nowroozi, 
1968). 

Sea  floor  noise  levels  at  Site  595  were  comparable  to  the  highest  values  yet  observed  by 
Scripps  OBS’s  in  the  Pacific  (Figure  14);  hence,  even  with  the  reduction  achieved  by  placing  the 
sensor  down  a  borehole,  the  MSS  noise  levels  are  on  the  order  of  30  dB  greater  that  those 
obtained  at  quiet  continental  sites  (Figure  15).  If  the  relationship  between  sea  floor  and  sub¬ 
bottom  noise  spectra  seen  in  Figure  14  is  maintained  at  quieter  locations,  however,  then  the 
down-hole  noise  levels  may  be  reduced  substantially  and  could  compare  favorably  with  typical 
continental  sites. 

Future  deployments  of  the  MSS  using  a  triaxial  set  of  seismometers  and  a  long-term  ocean 
bottom  recording  package  are  planned.  Such  a  system  complimented  by  an  array  of  OBS’s  should 
be  able  to  definitively  answer  questions  concering  the  mode  and  generation  mechanism  of  micro- 
seismic  noise  in  the  deep  ocean. 
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Figure  1.  Location  of  the  seismic  instrumentation  used  in  the  noise  experiment.  The  borehole 
instrumentation  package  of  the  Marine  Seismic  System  (MSS)  rested  undamped  at  the  bot¬ 
tom  of  Hole  395  A  (22'  45.35'  N  ;  46*  01.91  ’  \V),  drilled  five  years  earlier  on  Leg  45.  The 
ocean  bottom  seismometer  (OHS),  one  of  four  deployed  during  Leg  78H.  was  approximately 
1-2  km  west  of  Hole  395A.  The  OHS  launch  point  is  shown  (22  "  45.50'  N  ;  40  ’  05.43  '  W). 


Figure  2.  Schematic  depiction  of  the  MSS-data  channel  gains.  The  outputs  of  the  two  vertical- 
component  seismometers,  one  configured  50  cm  below  the  other,  were  split  into  3  channels 
each,  and  the  gains  of  the  resultant  six  channels  were  staggered  to  yield  144  dB  of  dynamic 
range.  lU’  and  lL’  refer  to  the  upper  and  lower  seismometers,  respectively.  The  use  of  a 
10-bit,  unsigned  digital  word  resulted  in  a  range  of  CO  dB  for  each  channel.  Adjacent  chan¬ 
nels  of  a  seismometer  overlapped  by  30  dB  and  upper  seismometer  gains  were  offset  from 
lower  seismometer  gains  by  24  dB. 
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Figure  3.  Displacement  responses  of  Leg  78B  seismic  systems.  The  MSS  employed  two  vertical- 
component,  short-period,  force-feed  back,  piezoelectric  seismometers  (Teledyne  (Jeotech 
model  S-700  ),  whose  outputs  were  filtered  to  prevent  aliasing  and  then  digitized  at  75  sam¬ 
ples  per  second.  The  MSS  response  is  peaked  near  11  Hz,  and  falls  at  12  dH/octave  towards 
lower  frequencies  and  at  15  dH/octave  to  the  Nyquist  frequency.  The  University  of  Texas 
(UT)  ocean  bottom  seismometer  system  consisted  of  a  triaxial  set  of  critically-damped,  10- 
Hz,  Mark  L6  gcophones  whose  outputs  were  sampled  with  an  8.352  millisecond  digitizing 
interval.  The  UT  OHS  response  increases  with  frequency  at  18  dH/octave  up  to  10  Hz,  and 
at  6  dH/octave  thereafter. 
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Figure  4.  Times  of  important  operations  and  events  during  the  Leg  78B  noise  experiment.  The 
borehole  instrument  package  (BIP)  First  entered  Hole  395A  at  0257Z,  28  March,  81,  but  was 
not  in  place  for  the  refraction  and  noise  experiments  until  0000Z,  29  March.  The  BIP  was 
recovered  at  0512Z,  30  March.  The  BIP  was  not  clamped  in  the  hole  but  was  isolated  from 
ship  motions  by  0.7  km  of  electro-mechanical  (E-M)  cable  payed  out  on  the  sea  floor  at 
1100Z,  29  March.  Twenty-six  hours  of  data  were  recorded  by  the  MSS  continuously, 
except  for  a  35-minute  gap  beginning  at  0545Z  29  March.  The  refraction  experiment  com¬ 
menced  at  0925Z,  29  March,  and  ended  at  0220Z,  30  March,  and  was  broken  by  a  two-hour 
hiatus  between  1800  and  2000Z,  29  March,  during  which  the  noise  experiment  was  per¬ 
formed.  The  University  of  Texas  ocean  bottom  seismograph  (UT  OBS)  recorded  eleven  34- 
second  noise  samples,  six  during  the  noise  experiment.  The  remaining  five,  denoted  by  dia¬ 
monds,  correspond  to  windows  when  refraction  shots  failed  to  explode. 


figure  5.  Sample  from  Channel  Ll  illustrating  characteristics  of  MSS  noise  data.  Low-frequency 
(0.2-0. 3  Hz)  energy  dominates  the  time  scries.  Features  which  contribute  to  noise  arc  the 
dropped  bit  near  11.5  seconds;  the  frequent,  random,  unidirectional  transients,  such  as  the 
one  labeled  near  10.5  seconds;  and  the  truncated  low-frequency  peak  centered  at  6  seconds, 
the  latter  possibly  be  due  to  an  instrument  non-linearity.  Units  are  digital  counts,  and  posi¬ 
tive  amplitude  corresponds  to  upward  mass  movement. 


Figure  6.  Power  spectral  densities  of  MSS  data  compared  with  estimates  of  instrument  noise. 
The  data  (solid  line)  were  computed  from  eight  13.65  second  samples  of  Channel  Ll  with  no 
instrument  correction  applied.  The  electronics  noise  was  estimated  from  RMS  noise  vol¬ 
tages  of  1/2  octave  bandwidth  provided  by  Tcledyne  Geotech.  The  least-count  noise  level  is 
that  expected  from  uniformly-distributed  errors  owing  to  the  representation  of  a  continuous 
quantity  by  a  finite,  digital  word.  The  transient-noise  spectrum,  which  is  attributed  to  the 
frequent,  impulsive  features  identified  in  Figure  5,  was  computed  from  equation  (A24)  and 
defines  the  high-frequency  limit  of  reliable  microseismic  noise  estimates.  Also  shown  in  this 
figure  is  the  squared,  normalized  MSS  displacement  response. 
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igure  8.  Coherences  between  the  highest  usable  gain  channels  of  each  of  the  MSS  seismometers 
(LI  and  U2).  The  coherence  spectrum  was  computed  from  twenty-four  109.2  second  sam¬ 
ples  of  each  channel.  Energy  between  the  two  seismometers  is  highly  incoherent  below  0.16 
Mr,  where  response  roll-off  causes  instrument  noise  to  dominate. 
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re  9.  Uncorrected  l-'T  OHS  power  densities.  Spectra  were  computed  from  eleven  17.1-s  seg¬ 
ments  obtained  from  six  34.2-s  noise  samples  recorded  over  a  75  minute  period.  The  large 
amplitudes  between  10  and  25  Hz  are  possibly  ascribable  to  resonance  of  the  OHS  frame, 
and  the  flat  portions  above  4  Hz  and  below  0.18  1h  are  likely  imposed  by  instrument  noise. 
Spectra  between  0.18  and  4  Hz  are  thought  to  estimate  reliably  ambient  noise  levels  on  the 
sea  floor. 


density  Inm^/Hx) 


0759.  29  M»r  81 


1800.  29  M»r  81 
0300.  30  Mar  81 


Fraquancy  (Hz) 


Figure  10.  Displacement,  power  densities,  in  nanometers  squared  per  Hertz,  in  the  Hole  395A  at 
approximately  10  hour  intervals.  Kach  curve  was  computed  from  twenty-four  109.2-second 
samples  from  Channel  U2  The  spectral  shape  of  the  noise  was  stable  over  this  period  but 
amplitudes  increased  with  time.  J’ower  densities  fall  at  approximately  HO  dH/decade 
beyond  the  peak  at  0.2-0.,'J  Hz. 
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Figure  12.  A  comparison  of  ambient  seismic  noise  at  and  below  the  sea  floor  using  corrected  UT 
OBS  and  MSS  data  recording  during  the  period  1800-1900Z,  29  March,  81.  The  UT  OBS 
data  were  six  34.2-s  samples,  and  the  MSS  data  were  six  109.2-s  samples;  the  former  were 
broken  into  eleven  17.  J-s  segments,  and  the  latter,  into  eleven  of  54.6-s  duration. 
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F.gure  13.  Comparison  of  corrected  UT  OBS  and  MSS  noise  observations  recorded  during 
refraction-experiment  windows  when  shots  failed  to  explode.  The  upper  set  of  curves  were 
computed  from  UT  OBS  data  and  the  bottom  from  MSS  data.  The  UT  OBS  was  pro- 
Jjr*rn™"Ll°  rccord  a  n>«ximum  of  10  seconds  per  event,  but  in  practice  only  34.2  seconds  of 
data  (4096  samples)  were  available  for  analysis.  Eight  segments  of  512  time  points  are  used 
to  calculate  estimates  for  both  the  UT  OBS  and  MSS  data. 
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Figure  14.  Comparison  of  Leg  78B  noise  spectra  with  OBS  observations  in  the  Pacific.  The  UT 
CBS  data  arc  the  same  as  those  in  Figure  12,  while  the  MSS  data  arc  the  earliest  (and, 
therefore  lowest)  of  those  in  Figure  10.  The  two  dashed  curves  are  extrema  of  observations 
made  with  the  Scr.pps  Institution  of  Oceanography  OBSs.  The  ‘Noisiest  Site’  was  southeast 
of  Hawaii  and  the  data  were  fifty-four  8-second  samples;  the  ‘Quietest  Site’  was  on  the  East 
l  acilc  Ivisc  and  the  data  were  thirty-one  4-sccond  samples. 
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Figure  15.  Comparison  of  the  sub-bottom  noise  spectrum  at  Site  S95  vtilh  spectra  at  two  quiet 
continental  sites.  The  Queen  Creek,  Arizona,  curve  was  taken  from  Melton  (1977),  who 
extracted  it  from  Fix  (1972).  The  Lajitas,  Texas,  data  were  recorded  at  an  extremely  quiet 
surface  site  during  a  time  of  low  winds  (Herrin, 1982). 
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Appendix  I 

The  noise  equation  and  its  transformation  in  horizontal  spatial  coordinates  are  derived  here. 
The  result  for  transformed  time  variation  is  well-known  and  will  be  staled  without  proof. 

Al.l  The  Noise  Equation 

Denote  (z,y,t)  by  the  vector  (;  then  the  functional  dependence  of  the  field  variables  may 
be  written  as  /  (£)  and  a  (£,z ),  for  example.  Also,  the  differential  element  dxdydt  will  be 
represented  by  d 

Expand  C„.  suppressing  the  z -dependence: 


c..(o  =  <  iJW.G^o-*,)/  u,)iifj(2Gr((o+(t-(t)/  (e*)i> 

(All) 

=  f fd(>d(2 G(eo-f1)Gr(Co+ei-Cz)<  /  Ui)f  (€.)> 

(A1.2) 

Let  — 

fi  +  fs  Then 

C..(0  -  J (6) 

(A1.3) 

Let  = 

Then 

C..(0  =  JdtsCn  (fs)[ fd(4Gtf4)Gr(f4+(-(s)l 

(A1.4) 

-  f  </(*?// ((*) r(f-6) 

(A1.5) 

i.e. 

=  r(f)  *  c„  (o 

(A1.6) 

G„(x,y,z,t)  =  fff  dzsdyidtariz-z^y-yyzS-tjC,,  (zs,ys,ts) 

(A1.7) 

Now,  because  /  (z,y,t)  and  G(z,y,z,<)  are  axisymmetric  we  can  write  C //  =  C/j  (r  ,t) 
and  T  =  r(r,z,t)  where  r2  =  z*+y*.  Then  the  above  convolution  is  axisymmetric,  that  is 
C..  =  C„  (r  ,z  ,t ).  In  plane  polar  coordinates,  this  convolution  takes  the  form 


oo  2r  oo 

C.,(r,z,t)=  J  dt„  £dtt  £r0drer{\r-r0\  ,z,t-l0)Cff  (r0,t0)  (Al.8) 

where  rr,  =  rr0  cos^ . 
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A1.2  Transformation  of  the  Noise  Equation 

C».  (r  ,2  ,t )  can  be  expressed  as  the  convolution  of  two  functions  T(r  ,z  ,t )  and  C /  /  (r  ,t ), 
each  of  which  are  autocorrelations.  The  integral  expression  for  Cuu  simplifies  considerably  if  we 
transform  it  over  both  time  and  horizontal  distance. 

The  Fourier  transform  pair  h  (< )  and  H  (w)  is  defined  as 


H(u)  -  fdt  t iu>  h(t) 

-00 

00 

MO  =  fd HW 


(A1.9) 


(A1.10) 


The  convolution  of  /  (t )  and  g  (/ )  has  transform  F(w)G(w);  their  correlation  has  transform 
F  (u;  )C*(  or),  where  the  dagger  denotes  the  complex  conjugate. 

A  similar  result  holds  for  the  (zeroth  order)  Bessel  transform.  Define  the  transform  pair 
A(r)  and  H(k)  as 


00 

Jf(k)  =  frdrj.(kr)h(r) 
00 

A  (r  )  =  £kdkJ.(kr)d(k) 


(Al.ll) 


( A 1 . 1 2) 


The  convolution  h  (r  )  of  /  (r  )  and  g  [r  )  is: 


it  oo 


/>(r)  =  '[d<fi0  Jf.dr,  /  (I  r-r,  |  )g{r0) 


(A1.13) 


and  their  cross-correlation  «  (r  )  is: 


7w  oo 


(r  ):  £dt0  £r0dr„  f  (|  r+r„  |  )g  (r„  ) 


( A 1 . 14 ) 


The  transform  of  h  (r  )  can  be  written  as 


it  00  fit  00 

H{k)=  27  [r»ir'^r |rrfr«*k7(|r-re  |  )  (A1.15) 

where  the  integral  representation  J0{kr)  =  — —  fe,krdd  has  been  used,  with  k  defined  such 

2* 

that  it  =  |  k  |  and  k  r  =  Arcos^.  Convert  the  bracketed  expression  in  equation  (A1.15)  to 


Cartesian  coordinates,  and  make  the  substitution  R  =  t-tc 
Y  =-  y  -y0  .  Then 


2»r  oc  r 

H(k)  = 

oo  2r 

=  2ir  F  (k  )  Jr0dr0g(r0)  £e'kr‘d< 


=  2itF(k)G{k ) 


Thus,  H{h)  =  2ir  F  (k  )G  (k  ). 


Proceeding  as  above,  S  (k )  may  be  reduced  to 


oo  2r 

S(k)  =  2 ,F(k)  -±-  £rBdrcg(r0)  je'" 


But  J„  (kr  )  =  J0  (~kr  ),  so 


S  (Jfc )  =  2nF(k)G(k)  =  H(k ) 


A1.3  The  Transformed  Noise  Equation 
Using  the  results  above,  we  can  write 


C „(k,z,u)  =  T(k,z,u)Cn  (k,u) 


where 


(r(*,*,w)),v  =  <?j,(t,ar,w)<7y(t,*,w) 
As  a  specific  example,  T„  ( k  ,z  ,w)  =  G/(k  ,z  ,u)G,  (k  ,z  ,w),  so 

{C..  (*  ,z  ,w)]„  =  |  G,  (fc  ,z  ,w)  |  7Cff 


-  ,  ,  ",  *  w  "  W*  V"*  L'*  k"*  fc"* 
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Appendix  II 

In  this  appendix  we  model  the  transient  noise  observed  in  the  MSS  data  as  a  compound, 
stationary  Poisson  process  (Feller,  1968)  and  derive  an  expression  for  its  power  spectral  density 
function  on  the  time  interval  0  ^  t  <  T .  We  define  the  process 


Ss(t)  =  E  **(0 


(A2.1) 


where  {A*  }  is  a  sequence  of  identically  distributed,  independent,  stochastic  variables,  and  N  has  a 
stationary  Poisson  probability  distribution 


P(N=n) 


(A2.2) 


The  expected  value  of  N  is 


<N>  =  E  »  P{N~n) 


=  AT  , 


(A2.S) 


(A2.4) 


so  that  the  constant  A  may  be  interpreted  as  the  mean  rate  of  event  occurrence. 
The  {Xk  }  are  assumed  to  have  the  form 


AT*  (t )  =  akh  ( t -tk )  , 


(A2.5) 


where  h(t)  is  the  transient  shape,  and  {a*  }  and  {t*  }  are  the  transient  amplitude  and  occurrence 
times,  respectively.  N,{a*  },  and  {<*  }  are  assumed  to  be  mutually  independent  random  variables; 
in  particular,  the  {tk  }  are  assumed  to  be  uniformly  distributed  over  the  time  of  observation,  T . 

We  wish  to  find  the  power  spectral  density  function  P$  (/  )  of  S n  (t)  at  frequency  f  ■  Ps  >s 
the  Fourier  transform  of  the  autocovariance  function  Rss  (r)  of  S n  : 


00 

Ps(f)-2j  ltss(r)e-i3'"di 


(A2.6) 
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*«(*)=<  (5w(t)-<  SN{i)>)  (S„{HT)-<SN(t-,r)>)  >  (A2.7) 

The  factor  of  2  in  equation  (A2.C)  is  due  to  our  normalization  of  Ps  (f  )i  in  which  the  total  power 
per  Nyquist  bandwidth  is  equal  to  the  variance  of  the  time  signal  (the  energy  of  the  negative  fre¬ 
quencies  is  folded  into  the  positive). 

Note  that  since  S#  is  stationary,  <  S/^(t)>  is  independent  of  t,  and  equation  (A2.7)  is 
equivalent  to 

Rss  (r)  =  <  SN  (t)Sfii  {t  -K)>  —  <  SN  >  *  .  (A2.8) 

We  will  evaluate  the  two  terms  of  this  equation  separately.  The  mean  value  of  S h  is 


<  SN>  =  £  P(N  =  n)<  £**>  =  £  P(JV  =  n)n<  Xk>  (A2.9) 

n  =0  k=  1  n  =0 

Since  N,  {ak  },  and  {tk  }  are  mutually  independent, 


<  sk>  =  <  y  >  <  xk> 


(A2.10) 


where  for  each  Xk  , 


T 

<  Xk>  =  <akh(t-tk)>  =<a>  £h(t- 


tk  )dtk 


(A2.ll) 


where  we  have  used  the  assumption  that  the  {ak  }  are  identically  distributed  and  therefore  have 
the  same  mean  <  a>  .  We  assume  that  the  interval  of  observation  T  is  much  larger  than  the 
transient  duration,  so  that,  approximately, 


<  Xk  >  = 


=  <«>  h  <x 


(A2.12) 


and  thus,  using  equation  (A2.4), 


<  SN>  =  A<  a>  £  h[t)dt  <  S  > 


(A2.1S) 


a 


9 


We  now  consider  the  bivariate  moment  in  equation  A8: 


<S„(t)S„[t+T)>  =  £  P(N=n)<  £**(«)  £  A'm(t+r) 

n  -  0  1  m  =1 

The  mutual  independence  of  the  {At  }  implies 


<  £  A'*  (t )  2  ATm  (t  -r)>  =  n  <  Xk  (t  )Xk  (t  +r)>  +  n  (n  -1)<  X  > 

k  ~  1  TO  =  1 


The  bivariate  moment  of  A'*  at  lag  r  is 


<  (/  )Ak  ( t  4r)>  =  <  (ak  h{t-tk  ))(ak  h  (t -tk  +r))> 


T 

=  <a2>  ±  +r )dt  ' 


*  <«J>  -j*  *  MO 


where  ‘  *  ’  denotes  cross-correlation.  Collecting  these  results, 


<  SN(t)SN(t+r)>  =  <  N  >  <  a2>  -j  h  *  h  (r)  +  <  N  (N -1)>  <  X  > 


From  equation  A2  we  Find  that 


<  N(N-1)>  =  (AT)3 


Thus,  from  equations  A12  and  A 13, 


<  SN  { t  )SN  (t  +r)>  =  A<  q3>  h  *  h  [r)  +  <  S  >  7  , 


so,  finally, 


Rss  (r)  -  A<  a3>  h  *  h(r) 
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It  is  now  a  simple  matter  to  determine  /‘5(f)  using  equations  A7  and  A22. 
result  from  Fourier  analysis  (llraeewell,  19G5)  is  that 


OO 


H(f)\  2 


where  11(f)  is  the  Fourier  transform  of  h(t).  Hence, 


Ps(f)  =  2A<  a2>  |  H({)\  2 


A  well-known 


(A2.23) 


(A2.24) 


